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Abstract 

We investigate magnetic coupling between a monolayer of prototype single-molecule magnets 
Mni2 and a ferromagnetic Ni(lll) substrate through S, using density-functional theory (DFT) 
and a DFT+C/ method. Our DFT and DFT+{7 calculations show that the Mni2 molecules favor 
antiferromagnetic coupling to the Ni substrate, and that they possess magnetic moments deviated 
from the magnetic moments of isolated Mni2 molecules. We find that the magnetic easy axis of 
the Mni2 on Ni (whole system) is dictated by that of the Ni substrate. The antiferromagnetic 
coupling is, dominantly, caused by superexchange interactions between the magnetic moments of 
the Mn and the Ni substrate via the S, C, and O anions. Our findings can be observed from x-ray 
magnetic circular dichroism or scanning tunneling microscopy. 

PACS numbers: 75.30.Et, 75.50.Xx, 75.70.-1, 71.15.Mb 
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I. INTRODUCTION 



Recently, tailoring static and dynamic magnetic properties of individual atoms and 
molecules by local environmental factors or local probes has drawn a lot of attention. This 
aligns with an effort to develop efficient devices for magnetic storage, spin transfer, molec- 
ular spintronics, or quantum computation, based on nanometer-sized magnetic atoms or 
molecules. One promising candidate for such devices is a group of single-molecule magnets 
(SMMs) consisting of several transition metal ions interacting through ligands. An indi- 
vidual SMM has a large magnetic moment and a high magnetic anisotropy barrier caused 
by spin-orbit coupling (SOC) and Jahn- Teller distortionsii^. In order to build devices in- 
cluding magnetic molecules, stable monolayers of magnetic molecules must be formed on 
various substrates, and mechanisms of interactions between molecules and substrates should 
be understood. A great progress has been made for deposition of various SMMs, such as 
prototype SMM Mni2^"-, Mng^, Fe^^, CryNii^, and TbPc2^^, on metallic, semiconducting, 
or superconducting substrates^. Electronic and magnetic properties of SMMs on substrates 
were characterized using x-ray absorption and photoemission spectroscopy^^^i^"— , scanning 
tunneling microscopy (STM)^>i^, /3-detected nuclear magnetic resonance^, and x-ray mag- 
netic circular dichroism (XMCD)^ii^. On the theoretical front, the electronic structure and 
magnetic properties of a SMM Mni2 monolayer adsorbed on Au were studied using density- 
functional theory (DFT)^^. The theoretical calculations^^ revealed that the Mni2 molecules 
are weakly coupled to the Au substrate, and that charge and spin transfer occurs from the 
Au substrate to the Mni2, mainly, through linker molecules. 

Despite the great advances made in recent years, studies of SMMs adsorbed on substrates 
are currently limited to non-magnetic substrates. An understanding of magnetic coupling 
between SMMs and ferromagnetic (FM) substrates is, however, of a great importance for 
device applications and for investigation of properties of SMMs. For instance, the magnetic 
coupling can be used to control the reversals of the magnetic moments of SMMs or of the 
magnetization of FM substrates, depending on the coercive field strength of the SMMs 
relative to that of the FM substrates, similarly to exchange-bias systems- . Additionally, 
the magnetic coupling can influence electron transport properties through SMMs bridged 
between FM electrodes. 

In this work, we examine magnetic coupling between a monolayer of SMMs Mni2 and 
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a FM Ni(lll) substrate, using DFT and a DFT+U method^ii^. We consider a structure 
where SMMs Mni2 are adsorbed on a Ni substrate via two S atoms. Our DFT and DFT+U 
calculations show that the Mni2 molecules prefer antiferromagnetic (AFM) coupling to the 
Ni substrate, and that their magnetic moments differ from the magnetic moments of isolated 
Mni2 molecules. We find that the magnetic easy axis of the adsorbed Mni2 is determined 
by the magnetic easy axis of the Ni substrate which depends on a thickness of the sub- 
strate. We also clarify mechanisms of the magnetic coupling by taking into account direct 
exchange coupling such as hybridization and an indirect exchange coupling mechanism such 
as superexchange interactions^^. We find that superexchange interactions prevail over direct 
exchange in the magnetic coupling. Our results can be observed in experiments such as 
XMCD or STM. 

The outline of this paper is as follows. We introduce our computational methods in 
Sec. II. We systematically construct a Mni2 monolayer adsorbed on a Ni substrate via S, as 
well as present properties of an isolated Mni2 and a bare Ni slab in Sec. III. We discuss the 
sign, magnitude, and mechanisms of the magnetic coupling between the Mni2 and the Ni 
substrate, and the electronic and magnetic properties of the system of interest in Sec. IV. 
Finally, we make our conclusion in Sec.V. 



II. COMPUTATIONAL METHODS 

We use two DFT codes, SIESTA - and VASP^°, to compute the electronic structure and 
magnetic properties of a Mni2 monolayer on Ni(lll). SIESTA uses localized numerical atomic 
orbitals as basis sets, while VASP uses plane waves. In both SIESTA and VASP calculations, 
we use the Perdew-Burke-Ernzerhopf (PBE) generalized-gradient approximation (GGA)^ 
for exchange-correlation potential. We do not include SOC in our calculations unless spec- 
ified otherwise. First, we discuss set-ups and parameters for SIESTA calculations. We use 
TrouUier- Mart ins pseudopotentials^'^ with scalar relativistic terms and core corrections for 



all the elements except H. We construct corresponding basis sets as discussed in Refs. |24jj25 |. 
For Mn, we include 3p orbitals in valence states. For Ni, we use a default basis set of DZP. 
We use a mesh cutoff of 400 Ry. We carry out self-consistent calculations until a density 
matrix converges to within 3 x 10~^. Second, we discuss set-ups and parameters for VASP 
calculations. We use projector-augmented-wave (PAW) pseudopotentials^ for all the ele- 
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FIG. 1: (Color online) (a) Side view of the whole structure: a SMM Mni2 adsorbed on a FM 
Ni(lll) slab of six atomic layers, via two S atoms, (b) Part of the interface region. The vertical 
distances, zi and Z2, are listed in Table HI (c) Top view (in the xy plane) showing the twelve Mn 
ions labeled, the O bonded to the Mn(l), Mn(2), Mn(6), and Mn(lO) sites, and the C bonded to 



the O and S. Modified from Ref. 



29]. 



ments. For Mn, we consider 3p orbitals as a semicore. For Ni, valance states consist of 4s 
and 3d orbitals. We use a kinetic energy cutoff of 400 eV. We perform self-consistent calcu- 
lations until the total energy converges to within 1 x 10~^ eV. For DFT+U calculations^^, 
we use an on-site Coulomb repulsion U term of 4 eV"^ for Mn d orbitals only. The value 
of U is selected in order to reproduce photoemission spectra measured for a crystal of Mni2 
molecules^^. 



III. WHOLE STRUCTURE 

We consider a structure whose unit cell is shown in Fig. [I](a). In the unit cell, a single 
Mni2 molecule is adsorbed on a FM Ni(lll) slab of six atomic layers, via two S atoms. This 
structure is referred to as 'whole structure.' We separately optimize geometries of a bare 
Ni slab (without a Mni2 molecule) and an isolated Mni2 molecule, combining the optimized 
geometries to construct the whole structure. 
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Let us discuss the geometry and properties of a bare Ni slab. We find that the equihbrium 
lattice constant of bulk face-centered-cubic (fee) Ni equals 3.500 A(3.522 A) using SIESTA 
(VASP), which agrees with experiments^. With the lattice constant, a Ni(lll) slab of 1 x 1 x6 
atoms is constructed, and its interlayer separations are optimized with in-plane separations 
fixed, until the maximum force is less than 0.03 eV/A for SIESTA (0.02 eV/A for VASP). The 
optimized interlayer separations are listed in Table [Tll The spin moments of the six atomic 
layers are, from the topmost to the bottommost layers, found to be (0.74, 0.68, 0.64, 0.64, 
0.68, 0.74) /iB, respectively, using SIESTA, while they are equal to (0.68, 0.68, 0.65, 0.65, 
0.69, 0.68) /iB, respectively, using VASP. 

We briefly review the geometry and properties of an isolated Mni2 molecule. More details 
can be found in Ref.[14i]. A standard Mni2 molecule, [Mni20i2(CH3COO)i6(H20)4], has the 
total magnetic moment of 20 /xb in the ground state, where the eight outer Mn spins (each 
Mn^"*" having 4 fi^) align antiparallel to the four inner Mn spins (each Mn^"*" having 3 /xb). 
For a standard Mni2 molecule, both the lowest unoccupied molecular orbital (LUMO) and 
the highest occupied molecular orbital (HOMO) arise from the Mn d orbitals. To reduce 
computational cost, we simplify a standard SMM Mni2 into [Mni20i2(HCOO)i6], where 
the CH3 groups in the standard Mni2 are replaced by H'^^, and the four water molecules 
are removedi^. This simplification does not affect the magnetic core of the Mni2. We 
substitute S atoms for the two H atoms closest to the Ni surface [Fig. [It^b)], such that the 
Mni2 molecule bonds to the Ni surface via the S. The bond length between the C and S 
equals 1.88 A (Fig. Wljo), Table H]). This slightly modified form of Mni2 is referred to as 
a S-terminated Mni2. For a S-terminated Mni2 molecule, the magnetic moments of the 
Mn ions remain unchanged, but each S atom has the magnetic moment of 1 /ib aligned 
antiparallel to the net magnetic moment of the twelve Mn ions. As a result, a S-terminated 
Mni2 molecule has the total magnetic moment of 18 /iB in the ground state. The magnetic 
moments listed in Table IIIII are computed by placing a sphere around each atom with the 
radius given in Table IIVI (The radii for O and S used in this work differ from those in 
Ref. [l3|.) The atomically resolved magnetic moments (Table UTTl) are lower than the actual 
values because contributions from inter-atomic regions are not included. For a S-terminated 
Mni2, the LUMO and HOMO originate from the S p orbitals. 

We now discuss the geometry of a unit cell of the whole structure. We determine a unit 
cell in order for the following conditions to be satisfied: (i) one Mni2 molecule must fit the 
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Ni surface area, (ii) neighboring Mni2 molecules interact weakly, and (iii) the total number 
of basis sets per unit cell should be computationally feasible. A monoclinic unit cell of 
17.33 X 17.33 X 36.00 (17.43 x 17.43 x 32.00 A^) suffices to meet the above conditions for 
SIESTA (VASP). The unit cell [Fig. [T](a)] consists of a Mni2 molecule attached to a Ni slab of 
7x7x6 atoms via two S atoms and a vacuum layer of 13.72 A (9.69 A) for SIESTA (VASP). 
The two S atoms bond to hollow sites of the Ni(lll) surface and to the C atoms close to 
the Mn(l), Mn(2), Mn(6), and Mn(lO) sites [Figs. [Tl^b),(c)]. The whole structure contains 
the total of 382 atoms, equivalent to 3450 valence electrons or 5552 orbitals. We sample 
3x3x1 /c-points for our DFT calculations. Considering the large system size, full geometry 
relaxations are difficult to achieve. Instead, we carry out partial geometry relaxations using 
SIESTA, where the Mni2, S, and Ni surface layer relax with the five Ni layers fixed. We find 
that the partial geometry relaxations do not alter our findings discussed below. We thus, 
henceforth, consider the unrelaxed whole structure only. For the unrelaxed geometry, our 
calculations show that large forces (the z components only) act mainly on the two S atoms 
and the six Ni surface atoms bonded to the S. They are, respectively, —2.35, —2.35, 1.58, 
1.27, 1.39, 1.44, 1.54, and 1.68 eV/A, using SIESTA (-1.52, -1.52, 0.93, 0.70, 0.69, 0.73, 
0.74, and 0.86 eV/A, using VASP). 

IV. RESULTS AND DISCUSSION 

A. Magnetic coupling type and strength 

We consider two spin configurations in the whole structure: (i) a FM configuration where 
the magnetic moment of the Mni2 is parallel to the magnetization of the Ni slab, and (ii) an 
AFM configuration where the magnetic moment of the Mni2 is now reversed such that it is 
antiparallel to the magnetization of the Ni slab. Our DFT calculations show that the AFM 
configuration has a lower energy than the FM one, and that the energy difference between 
the two configurations equals 39.5 meV (51.6 meV), using SIESTA (VASP). The exchange- 
correlation potential within the PBE GGA does not fully treat self-interaction corrections 
of the localized d states in the Coulomb potential. As a result, residual self- interactions 
induce slightly more diffuse d electrons, which entails overestimated exchange coupling. To 
examine the effect of self-interaction corrections, we apply a DFT-ft/ methodic to the two 
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FIG. 2: (Color online) Spin-polarized DOS projected onto (a) the Mn d orbitals, and (b) the Ni 
surface d orbitals, all orbitals of the S-terminated Mni2, the p orbitals of the O, C, and S atoms, 
and the C p orbitals bonded to the S atoms. The DOS is calculated for the AFM configuration. 
The majority (minority) spin is denoted as red dashed (black solid) curves. Notice the different 
vertical scales in (b). 



spin configurations with f/=4 eV for the Mn d orbitals. Our DFT+f/ calculations using 
VASP show that the AFM configuration has still a lower energy than the FM, and that the 
energy difference equals 46.3 meV. This energy difference does not differ much from that 
using DFT (Table |V|) . so that we, henceforth, discuss results obtained using DFT, unless 
specified otherwise. 



B. Electronic properties of the w^hole structure 

We investigate the effect of the Ni slab on the electronic properties of the S-terminated 
Mni2 in the AFM configuration. Using a Methfessel-Paxton'^^^ smearing parameter cr=0.01 
eV in VASP, we calculate the density of states (DOS) projected onto the Mn d orbitals, 
the Ni d orbitals, and the p orbitals of the O, S, and C atoms. As shown in Fig. [2](a), 
the Mn d orbitals broaden due to interactions with the Ni slab. Large broadening occurs 
for the d orbitals of the Mn(6) and Mn(lO) sites near 0.5 eV below the Fermi level, E{. 
Notice that the Mn(6) and Mn(lO) sites are closest to the Ni surface [Fig. [11(c)]. The O 
p orbitals strongly hybridize with the Mn d orbitals (Fig. [2]). In the band structure of Ni 
[top panel of Fig. El^b)], the minority-spin d bands intersect Ef, while the majority-spin d 
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FIG. 3: (Color online) (a) Charge density and (b) magnetic moment density integrated over the 
xy plane, as a function of z for the AFM (red filled circles) and FM (black empty squares) config- 
urations, the bare Ni slab (blue thin solid curve), and the isolated S-terminated Mni2 (green thick 
solid curve). The z coordinates of the Ni surface layer, S, C, and O sites are marked. The Mn sites 

X n 

locate out of the range. The S atoms are located at z = 0.09 A, as in Ref.|14i|. 



bands are completely filled. The S p orbitals and the C p orbitals bonded to the S atoms are 
delocalized due to the interactions with the Ni slab [third and fourth panels of Fig. El^b)]. 
The S p orbitals are spin polarized, and strongly hybridize with the surface Ni d orbitals 
than the C p orbitals do. 

Let us now discuss the charge distribution of the whole structure as well as calculate 
charge transfer between the Ni slab and the S-terminated Mni2 for the FM and AFM 
configurations. We perform our calculations using VASP. As shown in Fig. [3t^a), the charge 
distribution of the bare Ni slab has a long tail which is deeply penetrated into the region 
beyond the S. The charge distribution of the whole structure coincides with that of the 
isolated Mni2 near the mid-distance between the S and the C [Fig. [3t^a)]. To compute the 
charge transfer, we integrate over the xy plane the charge density of the whole structure, 
Pwhoie{x , y , z) , and the charge density of the isolated S-terminated Mni2, PsMM,isoi{x,y, z). 
We then take a difference between the two integrated charge densities: Ap{z) = /^^(pwhoic — 
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FIG. 4: (Color online) The difference in charge density between the whole structure and the isolated 
S-terminated Mni2, Ap(z), as a function of z for the AFM (red filled circles) and FM (black empty 
diamonds) configurations. The z coordinates of the S, the C and O closest to the S, and the Mn, 
are marked. 

PSMM,isoi) dxdy. The charge-density difference, Ap{z) (Fig. H]), peaks at the hnker molecules, 
the S, into which the charge distribution of the Ni slab is penetrated [Fig.[3](a)]. We calculate 
the charge transfer by an integration of Ap{z) over z. We find that some amount of charge 
is transferred from the Ni slab to the Mni2 molecule. The magnitude of the charge transfer 
is sensitive to the lower bound of the integration because the exact boundary between the 
Ni slab and the Mni2 is hard to be determined in the whole structure. For instance, when 
we integrate Ap{z) from z = —0.94 A to z = 11.90 A, the charge transferred equals 13.77 
(13.78) electrons for the AFM (FM) configuration. Here z = —0.94 A represents the mid- 
distance between the Ni surface layer and the S, and the S atoms are located at z = 0.09 A, as 
in Ref.kj. 
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FIG. 5: (Color online) The difference in magnetic-moment density between the whole structure 
and the sum of the bare Ni slab and the isolated S-terminated Mni2, Apm{z), as a function of z 
for the AFM (red filled circles) and FM (black empty diamonds) configurations. The z coordinates 
of the six Ni atomic layers, the S, the C and O closest to the S, and the Mn, are marked. 

C. Magnetic properties of the whole structure 



We examine how the adsorption modifies the magnetic moments of the S-terminated 
Mni2 and of the Ni slab as a function of z for the FM and AFM configurations, using VASP. 
Similarly to the charge-transfer calculation, we integrate over the xy plane the spin density 
of the whole structure, p^hoie; t^e spin density of the isolated S-terminated Mni2, PsliMisob 
and the magnetization of the bare Ni slab, p^j ^are- We refer to these integrated quantities 
as one-dimensional magnetic-moment densities and magnetization depending on z. We 
consider the following difference in one-dimensional magnetic-moment density: Apm{z) = 
LyiPwhoie - (LyPm^re ± Ly pfMM,isoi)]dxdy , whcrc the positive (negative) sign in front of 
the third term is applied to the FM (AFM) configuration. As shown in Fig. the most 
prominent effects of the adsorption in the difference Apm{z) are (i) that the difference peaks 
at the S ions, and (ii) that the AFM configuration possesses a smaller peak amplitude in 
the difference than the FM configuration does. For the isolated Mni2, each S ion has the 
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magnetic moment of 1 /ib, while for the whole structure, the S ions lose the spin polarization 
to great extent. As a result, the peak in ^pm{z) occurs at the S ions. Notice that in the 
isolated Mni2, the magnetic moment of each S ion aligns antiparallel to the net magnetic 
moment of the twelve Mn ions [Fig. El^b)]. In the whole structure, however, the magnetic 
moment of each S ion aligns parallel to the magnetization of the Ni slab, independent of 
the sign of the magnetic coupling between the Mni2 and the Ni slab (Table IIII|) . This 
entails that the peak amplitude for the AFM configuration is smaller than that for the FM 
configuration. (Effects (i) and (ii) are further discussed in Sec.IV.D.) For the Ni slab, the 
magnetic moments of the top two layers in proximity to the Mni2 substantially change due 
to the interactions with the Mni2. As shown in Fig. great changes appear near the sites 
where the S ions are bonded. The spatial dependence of the magnetic moments induced in 
the Ni slab for the FM configuration is similar to that for the AFM illustrated in Fig. |6l 
The main features discussed in this section remain unchanged with the application of the 
GGA+U method. 

Let us now discuss a way to calculate the total magnetic moment induced in the Ni 
slab, AM^j*^, and the magnetic moment of the adsorbed Mni2, M^^y^, within the GGA and 
GGA+U formalisms, using VASP. The difference in magnetic moment between the whole 
structure and the bare Ni slab, [M^hoie — M^f^], has contributions from both /S.M^^ and 
M|^Sj^, so that [M^hoic - = M|mm + AM^f ]. We first compute [M^hoie - M^^] 

including contributions from inter- atomic regions. To find the value of AM^j'^, we integrate 
the difference in one- dimensional magnetic-moment density between the whole structure and 
the bare Ni slab, /^^(Pwhoie "Pm.bare)'^^'^!/ (Fig. [7j), from the z coordinate of the bottommost 
Ni layer (farthest from the surface layer, z = —12.10 A) to the mid-distance between the Ni 
surface layer and the S atoms {z = —0.94 A). Similarly to the charge-transfer calculation, 
the value of AM^^ depends on the upper bound in the integration, and the difference in 
AM^^ between the AFM and FM configurations is very small (Fig. [7]). Then using the 
calculated values of [M^hoic - M^f ^] and AM^^, we find the value of M|^"m- 

Within the GGA, the difference, [M^hoie - amounts to -19.29 (20.54) /ig for the 

AFM (FM) configuration. When we integrate the one-dimensional density difference from 
z = -12.10 to ;z = -0.94 A, we find that AM^f amounts to 0.24 /ig (0.23 ^b) for the 
AFM (FM) configuration (Fig. E]). Thus, M|^V equals -19.53 (20.31) /^b for the AFM 
(FM) configuration. Within the GGA+U formalism, the difference, [M^hoie — -^Nr^] reaches 
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FIG. 6: (Color online) Magnetic moments induced in (a) the topmost Ni surface layer and (b) 
the second Ni layer from the surface in units of /xb for the AFM configuration. Not all of the Ni 
sites per atomic layer are shown. The dashed circles in (a) represent the sites where the S ions are 
bonded. The sites marked by the triangles in (a) correspond to the dashed circles in (b) in the xy 
plane. 



-19.05 (20.46) /iB for the AFM (FM) configuration. Since the U term is added to the Mn 
d orbitals only, the total magnetic moment induced in the Ni slab would not change with 
the addition of the U term. Thus, M^^y^ equals -19.29 (20.23) fiB for the AFM (FM) 
configuration. Considering that the magnetic moment of an isolated S-terminated Mni2 
equals 18 fiB, induced magnetic moments on the adsorbed Mni2 are shown in Table IVll 

To summarize the above calculations, when a S-terminated Mni2 molecule is adsorbed on 
a FM metal such as Ni, both the GGA and GGA+U calculations suggest that the adsorbed 
Mni2 molecule has the magnetic moment deviated from 20 fis, and that the sign of its 
deviation from 20 /i^ depends on the sign of the magnetic coupling between the Mni2 and 



12 



Ni Ni S CO 




Vertical coordinate z (A) 

FIG. 7: (Color online) The difference in one-dimensional magnetic-moment density between the 
whole structure and the bare Ni slab, /^.^(p^hoie ~ /'Ni bare)^^^^' ^ function of z for the AFM 
(red filled circles) and FM (black empty squares) configurations. 

the FM substrate (Table rvTj) . This result differs from the case that a S-terminated Mni2 
molecule is adsorbed on a non-magnetic metal such as Au. The magnetic moment of a 
S-terminated Mni2 adsorbed on Au increases to 19.8 /^b^- Notice that an isolated standard 
Mni2 molecule has the magnetic moment of 20 /i^. 

We now investigate the effect of the adsorption on the magnetic anisotropy of the AFM 
configuration. We perform DFT calculations including SOC using VASP. Considering the 
large system size, we treat the SOC non-self-consistently in coUinear cases. Previous the- 
oretical calculations^ have shown that a non-self-consistent treatment of SOC produces a 
slightly lower magnetic anisotropy barrier (lower by 14% for an isolated standard Mni2) 
than an experimental value or a self-consistently calculated value with a correct magnetic 
easy axis^^. Our calculation reveals that a bare Ni(lll) slab of six atomic layers possesses 
the magnetic easy axis parallel to the surface (which agrees with experiments^), and that 
the magnetic anisotropy barrier of the slab (7x7x6 atoms) equals 66.9 meV. An isolated 
S-terminated Mni2 molecule has, however, the magnetic easy axis normal to the surface, and 
its magnetic anisotropy barrier is found to be 4.88 meV. For the whole structure, we consider 
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only two coUinear spin configurations: (i) both the magnetic moment of the Mni2 and the 
magnetization of the Ni slab are parallel to the Ni surface, and (ii) they are perpendicular 
to the surface. Our calculations show that the first collinear configuration provides a lower 
energy than the second configuration, which implies that the Ni slab prevalently contributes 
to the magnetic anisotropy of the whole structure. The energy difference between the two 
collinear configurations equals 55.3 meV. This energy difference is lower than the magnetic 
anisotropy barrier of the bare Ni slab because for the given molecular orientation the mag- 
netic easy axis of the isolated S-terminated Mni2 is perpendicular to that of the bare Ni 
slab. A more quantitative analysis requires a self-consistent treatment of SOC and inclusion 
of noncollinear spin configurations and dipolar interactions, which is beyond the scope of 
the current study. 

D. Mechanisms of magnetic coupling 

To understand the mechanisms of the magnetic coupling between the Mni2 and the Ni 
slab, we consider direct exchange coupling and superexchange coupling-^. In the whole 
structure, the shortest separation between the Mn ions and the Ni surface atoms (6.58 A) is 
too large for direct exchange coupling to be effective. Thus, the Ni d orbitals do not hybridize 
with the Mn d orbitals. As a result, superexchange and/or Ruderman-Kittel-Kasuya-Yosida 
(RKKY)-like interactions^!^ may cause the magnetic coupling. In numerous systems, either 
superexchange or RKKY-like interactions play a role in magnetic coupling. In this work, 
we consider superexchange coupling exclusively because it explains the magnetic coupling 
of interest. 

Let us examine contributions of superexchange interactions to the magnetic coupling. In 
the whole structure, the shortest separation between the S sites and the Ni surface atoms 
is set to be 2.52 A [Fig. [I](b), Tabled]. As discussed, charge transfer occurs from the Ni 
slab to the S-terminated Mni2 molecule, or to the S ions which is the closest to the Ni 
slab [Fig. El^a)]. The fact that the LUMO and HOMO of the Mni2 arise from the S p 
orbitals, corroborates the charge transfer to the S ions. The band structure of Ni entails 
favorable transfer of minority-spin electrons. Thus, the total magnetic moment induced in 
the Ni slab, AM^;'^, is positive (Sec.IV.C), and the difference in one-dimensional magnetic 
moment density between the whole structure and the bare Ni slab (Fig. [7j) has a positive 
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FIG. 8: Magnetic coupling via superexchange interactions for (a) the whole structure and (b) a 
single Mn atom adsorbed on Ni(lll) via S. Magnetic coupling via hybridization for (c) a single Mn 
atom directly bonded to Ni(lll). The majority and minority-spin Ni 3d bands at zero temperature, 
S 2p orbital levels (degenerate), and Mn 3d orbital levels, are illustrated. The Mn d levels for (b) 
and (c) are degenerate. 

peak amplitude at the Ni surface layer. This charge transfer lowers the magnetic moments 
of the S ions, resulting in the magnetic moments parallel to the magnetization of the Ni 
slab [Fig. [Ht^a)]. Let us now focus on the region between the S and the Mn ions [Fig. [11(b)] . 
The magnetic moments of the C and O ions closest to the S ions are much smaller than 
those of the S ions (Table HTB . As shown in Fig. [Hb), C and O anions are located between 
the S and the Mn(l) and Mn(6). The bond angle created by the Mn(6) [or Mn(l)] and 
the O and C bonded to the Mn [Fig. [^b)] equals 123 degrees. The bond angle of S-C-0 
in Fig. lU^b) falls on 114-118 degrees. While the S ions have dominant contributions from 
2px and 2py orbitals, the C and O anions have equal contributions from 2px, 2py, and 2pz 
orbitals. The Mn(l) (Mn^+) possesses exactly half-filled spin-polarized t2g orbitals. Our 
calculations suggest that this inner Mn ion is ferromagnetically coupled to the magnetic 
moments of the S, C, and O ions. The Mn(6) (Mn^+) has spin-polarized half-filled t2g-like 
orbitals and quarter-filled Cg-like orbitals. Our calculations indicate that the outer Mn ion is 
antiferromagnetically coupled to the magnetic moments of the S, C, and O ions. The same 
logic is applied to the other linker molecule. Therefore, the magnetic moments of the S, C, 
and O ions are antiferromagnetically coupled with the net magnetic moment of the twelve 
Mn ions. As a result, the superexchange interactions render overall AFM coupling between 
the local magnetic moments of the Ni and the total magnetic moment of the Mn ions via 
the S, C, and O p orbitals [Fig. |8]^a)]. We emphasize that not only the separation between 
the Ni and the Mn ions but also the type and geometry of bonding to the Ni and Mn, play 
an important role for the superexchange interactions. 
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Ni Ni S Mn 




Vertical coordinate z (A) 

FIG. 9: (Color online) The difference in one-dimensional magnetic-moment density between system 
(i) and the bare Ni slab, /^..^ (p^stem (i) ~ /'Ni barc)^^'^^' ^ function of z for the AFM (red filled 
circles) and FM (black empty squares) configurations. Calculated using VASP. 

To gain further insight into the nature of the magnetic coupling, we investigate magnetic 
couphng in the following testbed systems: (i) a single Mn atom adsorbed on a Ni(lll) slab 
via one S atom [Fig.[8](b)], and (ii) one Mn atom directly bonded to a Ni(lll) slab [Fig.[8](c)]. 
The magnetic coupling in system (i) seems to have similar mechanisms to that in the whole 
structure, in that the magnetic moment of the Mn interacts with the Ni slab through the 
S, and yet system (i) is much simpler than the whole structure. Now if we remove the S 
atom from system (i) and bring the Mn atom closer to the Ni slab, the characteristics of the 
magnetic coupling are qualitatively modified. That is the case for system (ii). 

Let us first discuss the magnetic coupling in system (i). We set the shortest distance 
between the S and Ni surface atoms to be the same as that for the whole structure (Table [T]), 
and assign positive spin polarization to the Ni slab in both AFM and FM configurations. 
The S ion bonds to hollow sites of a Ni slab of 5 x 5 x 6 atoms, and the bond angle of 
Mn-S-Ni equals 145.5°. Our DFT calculations using SIESTA suggest that similarly to the 
whole structure, the Mn ion favors AFM coupling to the Ni slab via the S anion (Table [11). 
As in the whole structure, a major contribution to the magnetic coupling originates from 
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a superexchange interaction via the S anion. An isolated Mn-S dimer has the ground-state 
spin of S = 5/2, where the Mn ion carries most of the spin polarization. Since the S p 
orbitals are completely filled and the minority-spin Ni d bands are not fully filled, minority- 
spin electrons can be transferred now from the S to the Ni slab [Fig. [H](b)]. Consequently, 
the S p orbitals possess slight spin polarization parallel to the magnetization of the Ni slab, 
and negative magnetic moments are induced in the Ni slab (Fig. [9]). The S p orbitals are 
antiferromagnetically coupled to the Mn d orbitals. Overall, this leads to AFM coupling 
between the Mn and the Ni (Fig. Mjo), Table H]). The magnetic coupling in system (i) is 
stronger than that in the whole structure due to the following reasons: (1) The vertical 
separation between the Mn and the Ni surface atoms in system (i) is shorter than that in 
the whole structure (Table H]). (2) A superexchange interaction through S is stronger than 
that via O or C, since S atoms more diffuse. In the whole structure, the superexchange 
interactions are weakened by pathways through the S, C, and O orbitals. For the given 
geometry of system (i), whether the 3d orbital levels of the adatom are less than or more 
than half-filled, would not change the sign of the magnetic coupling because the magnetic 
coupling is mainly governed by a superexchange interaction. 

We now examine the magnetic coupling in system (ii). In this system, the Mn atom 
directly bonds to hollow sites of a Ni slab of 5 x 5 x 6 atoms, and the Mn d orbital levels 
are exactly half-filled [Fig. Et^c)]. Our DFT calculations using SIESTA show that stable FM 
coupling is formed between the Mn and the Ni (Table H]), in contrast to system (i) and 
the whole structure. In system (ii), the bond length between the Mn and the Ni is short 
(Table [Tj), and so the minority-spin Mn d orbitals hybridize with the minority-spin Ni d 
orbitals. This results in FM coupling between the Mn and the Ni, which agrees with other 
theoretical calculations^^. In system (ii), if the 3d orbital levels of the adatom are less than 
half-filled, then hybridization favors AFM coupling between the adatom and the Ni^. 

V. CONCLUSION 

We have investigated magnetic coupling between the SMM Mni2 and a FM Ni(lll) 
substrate, using DFT and a DFT+U method. Our DFT and DFT+U calculations have 
shown that the Mni2 interacts antiferromagnetically with the Ni substrate via S, and that 
the magnetic moment of the Mni2 is deviated from that of an isolated standard Mni2, due to 
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interactions with the Ni substrate. We have also found that the magnetic anisotropy of the 
whole structure is dictated by the anisotropy of the Ni slab, and that the magnetic anisotropy 
barrier of the whole structure depends on the orientation of the Mni2 relative to the Ni slab. 
The magnetic coupling is mainly caused by superexchange interactions between the magnetic 
moments of the Ni and Mn via the S, C, and O anions. Superexchange interactions take 
into account the bonding types and bonding geometries of linker molecules to the Ni slab 
and the Mni2. Our findings can be observed in XMCD experiments and provide insight into 
magnetic couphng of other types of SMM to FM substrates. 
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TABLE I: Vertical distances z (in A), bond lengths d (in A), and an energy difference Ai? (in eV) 
between the AFM and FM configurations. < (A£' > 0) implies stable AFM (FM) coupling. 
For the calculations of IS.E using SIESTA, we sample 3x3x1 fc-points for the whole structure and 
4x4x1 fc-points for systems (i) and (ii). 





zx = z(Ni-S) 


Z2 = z(S-Mn) 


d(Ni-Mn) 


d(Ni-S) 


d(S-C) 


A^ 


whole structure 


2.08 


4.48 


6.58 


2.52 


1.88 


-0.0395 


system (i) 


2.08 


2.40 


4.70 


2.52 


N/A 


-0.1050 


system (ii) 


N/A 


N/A 


2.80 


N/A 


N/A 


0.1692 



TABLE II: Optimized interlaycr separations (in A) for the Ni slab of six atomic layers. Azjj 
indicates the vertical distance between the ith and the jth layers, a is the equilibrium lattice 
constant for bulk fee Ni. 





SIESTA 


VASP 


AZ12 


2.011 


2.008 


AZ23 


2.028 


2.037 


AZ34 


2.021 


2.039 


AZ45 


2.023 


2.029 


AZ56 


2.011 


2.007 


a 


3.500 


3.522 
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TABLE III: Atomically resolved magnetic moments for an isolated S-terminated Mni2 and the 
AFM and FM configurations, in units of Bohr magneton (/Us), within the GGA, using VASP. The 
numbers in the parentheses denote magnetic-moment differences between the isolated Mni2 and 
the whole structure, where only the magnitudes are listed. The Mn ions are labeled in Fig. HJc). In 
the last row, the total magnetic moment of the Mni2 is calculated from the sum of the atomically 
resolved magnetic moments. 



AtoTTiir sDPrips 




whole AFM 


whole FM 


iviii 1 1 j oi ivin I ^ ) 


9 R'i 
— i^.UO 


9 fi9 (0 ni 


9 fin (C\ D'^'l 


Mnl'.S') or MnM") 


—2.59 


2 59 (() on) 


—2 58 ("0 Oil 


Mn(5) or Mnl'9') 




— 54 (T) ni 


3 54 fn 01 


Mn(7) or Mn(ll) 


3.52 


-3.54 (0.02) 


3.54 (0.02) 


Mn(6) or Mn(lO) 


3.50 


-3.56 (0.06) 


3.57 (0.07) 


Mn(8) or Mn(12) 


3.56 


-3.57 (0.01) 


3.57 (0.01) 


Mn (total) 


17.78 


-18.00 (0.22) 


18.05 (0.27) 


closest to S 


-0.01,-0.01,-0.04,-0.04 


0.00, 0.00, 0.00, 0.00 


0.01, 0.01, 0.02, 0.02 


(total) 


0.33 


-0.44 


0.53 


C closest to S 


0.01, 0.01 


0.01, 0.01 


0.00, 0.00 


C (total) 


0.21 


-0.18 


0.19 


S(l) or S(2) 


-0.44 


0.09 


0.10 


Mni2 (total) 


17.44 


-18.41 


18.95 
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TABLE IV: The atomic sphere radii used for calculations of atomically resolved magnetic moments 



using VASP (Table imi). 



Atomic species 


Atomic radius (A) 


Ni 


1.29 


Mn 


1.32 


S 


1.16 





0.82 


c 


0.86 


H 


0.37 



TABLE V: Energy difference IS.E (in eV) between the AFM and EM configurations, within DET 
and DET+[/ formalisms, using SIESTA and VASP. 



GGA, SIESTA 


GGA, VASP 


GGA+;7, VASP 


-0.0395 


-0.0516 


-0.0463 



TABLE VL Total and induced magnetic moments (in ^ub) of the adsorbed Mni2 for the AEM and 
EM configurations, within the GGA and GGA+U formalisms, using VASP. The induced moments 
are obtained from comparison with the magnetic moment of an isolated S-terminated Mni2. 





Total (GGA) 


Induced (GGA) 


Total (GGA+U) 


Induced (GGA+U) 


AFM 


-19.53 


-1.53 


-19.29 


-1.29 


EM 


20.31 


2.31 


20.23 


2.23 
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